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ABSTRACT 

The final masses of Jovian planets are attained when the tidal torques that they exert 
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o 
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, on their surrounding protostellar disks are sufficient to open gaps in the face of disk 

bJOl viscosity, thereby shutting off any further accretion. In sufficiently well-ionized disks, the 



predominant form of disk viscosity originates from the Magneto-Rotational Instability 
(MRI) that drives hydromagnetic disk turbulence. In the region of sufficiently low 
ionization rate - the so-called dead zone - turbulence is damped and we show that lower 
mass planets will be formed. We considered three ionization sources (X-rays, cosmic 
, rays, and radioactive elements) and determined the size of a dead zone for the total 

^f) i ionization rate by using a radiative, hydrostatic equilibrium disk model developed by 

I Chiang et al. (2001). We studied a range of surface mass density (So = 10^-10^ g cm I 

' and X-ray energy [kTx = 1 — 10 keV). We also compared the ionization rate of such a 

■ disk by X-rays with cosmic rays and find that the latter dominate X-rays in ionizing 
protostellar disks unless the X-ray energy is very high (5 — 10 keV). Among our major 

Q_i' conclusions are that for typical conditions, dead zones encompass a region extending out 

Q . to several AU - the region in which terrestrial planets are found in our solar system. Our 

results suggest that the division between low and high mass planets in exosolar planetary 
^ I systems is a consequence of the presence of a dead zone in their natal protoplanetary 

disks. We also find that the extent of a dead zone is mainly dependent on the disk's 

■ surface mass density. Our results provide further support for the idea that Jovian 
^ • planets in exosolar systems must have migrated substantially inwards from their points 
^ of origin. 

Subject headings: circumstellar matter - planetary system - solar system: formation - 
accretion disks - MHD - stars: pre-main-sequence 



1. Introduction 

The discovery of exosolar Jovian-mass planets has stimulated intense efforts to understand 
how planets form (e.g. Lin et al. 1999; Wuchterl et al. 2000). The current theoretical models 
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are based on two general types of physical processes - disk instability or core accretion. In the 
disk instability picture (e.g. Cameron 1978; Boss 1997; Mayer et al. 2002), gravitational instability 
in a sufficiently massive disk fragments it into the clumps out of which planets may form. This 
process may require as little as ~ 10^ years to form Jovian planets, but may not be able to explain 
sub- Jupiter mass planets (Boss 2000). This is in stark contrast to the core accretion model (e.g. 
Mizuno 1980; Pollack et al. 1996) wherein the coagulation of dust grains ultimately leads to the 
formation of planetesimals that collide with one another to form planetary cores. When a core 
reaches a critical mass in the range ~ 10 — 15M0, gas around the protoplanet quickly accretes onto 
it ("runaway accretion") forming a gas giant. This process takes at least ~ 10^ — lO'^ years - the 
average disk lifetime (e.g., Feigelson &; Montmerle 1999), and may not be able to explain massive 
Jovian planets. 

Our paper focuses on an important aspect of Jovian planet formation that is common to both 
pictures. Regardless of whether proto- Jovian planets form through core accretion or gravitational 
instability, accretion onto a planet will continue until it is massive enough for the tidal torque that 
it exerts upon the surrounding disk to overcome the viscous disk torque which continues to fill-in 
any gap (e.g. Lin & Papaloizou 1985). Planets that are sufficiently massive to open gaps in the 
disks have final masses that depend upon the disk viscosity (Nelson et al. 2000). In absence of 
viscosity, the tidal torque induced by the protoplanet opens a gap sooner thereby terminating the 
accretion earlier and leaving the planet with a smaller mass. 

Thus the final mass of Jovian planets comes down to the question of the origin of the viscosity 
of protostellar disks at the time that significant protoplanetary cores are present within them. If 
disks are turbulent, then turbulent viscosity is likely to be the major factor in determining these 
masses. The Magneto-Rotational Instability (MRI) is the most promising source of turbulent disk 
viscosity (Balbus & Hawley 1991). Even in the presence of an initially very weak magnetic field, 
the MRI will generate a significantly magnetized disk. This MRI viscosity requires good coupling 
between the gas and the magnetic field. Poor coupling - which occurs when the gas is insufficiently 
ionized - damps out the MRI and leads to a so-called "dead zone" where the viscosity is effectively 
zero (Gammie 1996). The high column density for protoplanetary disks ensures that dead zones 
are expected to occur during planet formation (Sano et al. 2000). 

Disk viscosity will not entirely vanish, even in the dead zone, however. The gravitational 
interaction between the fairly massive protoplanetary cores - the forerunners of the Jovian planets 
- and the disk generates density waves within the disk which can transport angular momentum and 
provide an effective "viscosity" in the disk. This disk-planet interaction leads to both planetary 
migration (Goldreich &; Tremaine 1980) as well as evolution of the disk through the nonlinear 
damping and shocking of the waves (Larson 1989; Goodman & Rafikov 2001, - GROl). The 
deposition of the angular momentum back into the disk as a consequence of this shock-induced 
wave damping is equivalent to the action of a viscous process, (Spruit 1987; Larson 1989) with an 
equivalent dimensionless a parameter of the order ct^amp — 

10-^-10-3 (e.g. GROl). This viscosity 
parameter is one to two orders of magnitude smaller than is expected for MRI turbulence in the 
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regions of the disk where the magnetic field is well coupled to the fluid. 

We can readily compare the masses of planets that are predicted for each of these mechanisms of 
disk viscosity. For the case of a turbulent viscous disk (e.g. Bryden et al. 1999; Rafikov 2002), the gap 
will open when the angular momentum transfer rate by tidal torque Ht = 0.23{Mp/M^yTiapQ'^ {ap/ h) 
exceeds that driven by viscous torque Hi, = SnTiuapO,. Here is the Keplerian angular velocity 
and v is the viscosity. Expressing the viscosity in terms of a standard a. parameter v ~ aturh^sh 
(Shakura &; Sunyaev 1973), the ratio of a gap opening planetary mass Mp to a stellar mass M* is 
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where the disk height and the radius of the planet's orbit are h and Op, respectively. 

In the case that disk viscosity is provided solely by damped density waves, then the ratio of 
the planetary masses has been calculated by Rafikov (2002); 
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where the Toomre parameter \s Q = i^Cs / {ttGT,) , Cs is the sound speed, and S is the surface mass 
density. Note that the gap opening condition for an inviscid disk by Lin &; Papaloizou (1985) 
corresponds to Q = 1 case. For the standard disk model by Chiang et al. (2001), we calculate that 
Q = 54 and h/up ~ 0.04 at 1 AU, so that the planetary mass is very small. 

By taking the ratio of these two predicted planetary mass scales, we immediately see that the 
planetary masses in well coupled zones in which the MRI is active will be much larger than in the 
dead zone for which only the oidamp value pertains. As an example, at 1 AU, and for a Toomre 
Q value of 54 calculated above, the mass of a Jovian planet in a well-coupled zone dictated by 
MRI turbulence Mp^t^rb will exceed the mass that a Jovian protoplanet will attain in a dead zone 
through density wave damping Mp^^amp provided that aturh exceeds a minimum value of; 



Oiturh > Oimin = min 
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which is less than lO^'^ for the aspect ratio at 1 AU {h/ap ~ 0.04). Numerical simulations of the 
MRI instability show that is much larger than this - by two orders of magnitude or more. 

Thus, a distinct jump in planetary masses - by a factor that is similar to the ratio of the masses 
of the Jovian to terrestrial planets - is expected at the disk radius that determines the maximum 
radial extent of the dead zone. 

In this paper, we perform a detailed investigation in order to calculate precisely where the dead 
zone occurs in well constrained protostellar disk models. This allows us to determine where in the 
disk that this "jump" in planetary masses is likely to occur. In our second paper (to be submitted), 
we calculate precisely what these planetary masses are. We proceed by carefully computing the 
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ionization balance within the protostellar disk, since it is this that determines the couphng of the 
magnetic field to the disk and hence, where the MRI induced turbulence can be sustained. The 
main sources of ionization in a circumstellar disk arc X-rays from the stellar plasma (probably 
produced by magnetospheric reconnection events) as well as cosmic rays from the outside of the 
stellar system. Glassgold et al. (2000) showed that X-ray ionization dominates cosmic ray ionization 
out to ~ 1000 AU in an optically thin disk. We compared these ionization sources for the self- 
consistent disk model developed by Chiang ct al. (2001). We include radioactive elements as a 
global and constant ionization source. Wc used the disk model developed by Chiang et al. (2001) 
primarily because it reproduces the spectral energy distributions (SEDs) of two T-Tauri and three 
Herbig Ae stars extremely well. 

Our major finding is that dead zones typically encompass the region where terrestrial planets 
are found in our solar system (< 2.0 — 3.3 AU). The physical definition of a dead zone is reviewed in 
§2. We introduce the disk model of Chiang et al. (2001) and elucidate the connection between the 
disk model and the X-ray ionization rate (using Glassgold et al. 1997) in §3. We state our results 
in §4 and apply them to the well studied case of AA Tau in §5. The important symbols used in 
this paper are summarized in Table 1. 



2. Dead Zones in Disk Models 

A dead zone is the region in a disk where MRI turbulence is unsustainable due to poor coupling 
of the field to the disk. The absence of turbulence is what gives this region its negligible viscosity. 
Gammie (1996) showed that MRI turbulence will be damped at a physical scale of A whenever the 
local growth rate of the MRI instability (~ Va/A) is balanced by the Ohmic diffusion at that scale 
(2± 77/A^). For complete damping at some radius of the disk, one requires that all the turbulence 
scales less than the local pressure scale height h{a) are damped; i.e. A < /i. By equating these two 
time-scales, one finds that the MRI turbulence damps in regimes with a small magnetic Reynolds 
number: 

ReM = ^<l. (4) 

ri 

Here, Va = B/^Attp ~ a\!^^Cs is the Alfven speed, p is the mass density, Cg = \/ kT/ /jluih is the 
sound speed {k is the Boltzmann constant, T is the temperature, p is the mean molecular weight, 
and niH is the mass of the hydrogen) and aturb is the viscosity parameter that measures the strength 
of the turbulence (Shakura & Sunyaev 1973). This formula is one of the most important in the 
paper because it is the condition that determines the dead zone region in a disk once the disk 
structure and ionization state are specified. The importance of the latter is readily established by 
noting that the Reynolds number in this equation is sensitive to the ionization of the disk, through 
its dependence on the diffusivity of the magnetic field, rj, which takes the form; 

Tj = — cm2 s-1 , (5) 
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Table 1: Some Important Symbols 
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Ionization rate 
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X-ray ionization rate 
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Cosmic ray ionization rate 
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Radioactive elements ionization rate 
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Total photoelectric absorption cross section 
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Optical depth 
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Distance from the X-ray source to the disk 
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Number density of neutral atoms 
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Nh 


Surface number density along the path of X-rays 


eq. 


(13), (15) 
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Vertical surface number density: A^_l = n{a, z')dz' 
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(15) 


X 


Vertical surface mass density: x = Pio-i z')dz' 


eq. 
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a 


Grazing angle from the central star 
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(21) 


a' 


Grazing angle from the X-ray source at (12i?0, 12i?0) 


eq. 


(15), (21) 
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Table 2: Parameters used in the models. 
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Table 3: Dead zone radii for the star with the X-ray luminosity of 10 erg s . The range of dead 
zone radius corresponds to the range of surface mass density. 
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where T is the disk temperature that is obtained from disk models (either interior Tj or disk surface 
Tils depending on the disk height we are interested in - see §2.3, eq. (3) and (4) in Chiang et al. 
2001). The electron fraction is = rie/n where rig and n arc the number density of electrons and 
neutral atoms respectively. The lower the electron fraction, the higher is the diffusivity coefficient 
for the field and hence the smaller is the magnetic Reynolds. Thus, poorly ionized disks will tend 
to have lower values of ReM, resulting in larger dead zones. 

One of the recent numerical MRI disk simulations by Fleming et al. (2000) defined the magnetic 
Reynolds number as 

R'eM ^ ' 

and determined the critical value R'^m crit ~ assuming the presence of a uniform vertical 
magnetic field. Another recent analysis by Balbus & Terquem (2001) suggested that the Hall effect 
should be included into the diffusion equation. The Hall effect is due to the velocity difference 
between electrons and ions, and causes a small transverse potential difference. The effect also 
diffuses the magnetic field just as the ohmic diffusion does. It implies that the magnetic Reynolds 
number might be smaller than the case of Fleming et al. (2000). 

Summarizing these results, it appears that the critical value of magnetic Reynolds number 
ReM,crit IS somcwhat uncertain and that there are a few different definitions for ReM- Taking these 
uncertainties into account, we used the critical magnetic Reynolds number of 1 and 100 with the 
alpha parameter of 1, 0.1, and 0.01. The extreme cases of {ReM,critiOtturb) =(1) 1)) (1) 0.01), (100, 
1), and (100, 0.01) correspond to R!^j^^ = 1, 10, 100, and 1000 respectively in the definition by 
Fleming et al. (2000). The following is the criterion we used to determine the dead zone: 

ReM = ^ ReM,crit ■ (7) 

The dead zone is the region in the disk where the magnetic Reynolds number RgM is smaller than 
some critical value ReM,crit = 1 or 100. This region is characterized as having no MRI turbulence. 

The calculation of the needed electron fraction Xg (see eq. (5)) is well established, at least in 
principle. Assuming steady state ionization balance, Oppenheimer &; Dalgarno (1974) derived an 
equation for the electron fraction in the gas; namely, 

(3d Pan (3d/3rn 

where xm is the metal fraction. In our case, ^ represents the ionization rate of X-rays, cosmic 
rays or radioactive elements, or total ionization rate of these three sources. Their mathemati- 
cal forms are found in §3. Three beta terms represent three kinds of recombination rate coef- 
ficients - the dissociative recombination rate coefficient for electrons with molecular ions {Pd = 
2 X 10~^T~^/^ cm^ s~^), the radiative recombination coefficient for electrons with metal ions 
{Pr = 3 X 10~^^r~^/^ cm^ s~^) and the rate coefficient of charge transfer from molecular ions 
to metal atoms (/3t = 3 x 10~^ cm^ s~^). 
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In this paper, we only take account of the two extreme cases - the metal poor (xm = 0) and 
the metal dominant {xm S> Xg). In both cases, the electron fraction is written as 




where the number density n depends on disk models (see eq. (17)) and /3 is a recombination rate 
coefficient. In the metal poor case, only the recombination of electrons with molecular ions (disso- 
ciative recombination) is important and /? = /J^ = 2 x 10~^T~^/^ cm^ s"-*^. We use this coefficient 
throughout the paper except §5. In the metal dominant case, only the recombination of electrons 
with heavy metal ions (radiative recombination) becomes important and the recombination rate 
coefficient is replaced by the radiative recombination coefficient /3 = /?,. = 3 x IQ^^^T^^^"^ cm^ s~^ 
(Fromang et al. 2002). The astrophysical significance of these two cases is related to the density 
regime of the protostellar disk. We compare the results for these two cases as applied to AA Tau, 
in S5. 



3. Ionization Rate of a Self-Consistent Disk Model 

The previous section outlined precisely how the dead zone within a disk can be calculated once 
a disk model and ionization rate can be specified. In this section, we review the three sources of 
ionization for protostellar disks - X-rays from the central protostar, cosmic rays, and radioactive 
elements that are mixed with the disk material - and then calculate their respective ionization rates 
in the context of our chosen disk model. There is one caveat concerning cosmic ray ionization of 
protostellar disks - cosmic rays are likely to be swept away from an region containing a turbulent 
MHD jet or outflow (e.g. Skilling & Strong 1976; Cesarsky Sz Volk 1978). Since jets are manifest 
in Class I and II Young Stellar Objects (YSOs), it is possible that the bulk of the ionization of 
these sources is achieved solely by their YSO X-ray fluxes. For completeness, we show both sets of 
results - ionization from purely X-ray as well as combined ionization. 

The combined contribution of each of these ionization sources produces the total ionization 
rate which is their sum; 

Ctot = Cx + CcR + Cra ■ (10) 

We summarize the calculation of each of these three rates below and then also compare the effects 
of X-rays with that of cosmic rays in our disk model calculations. 



3.1. Ionization Rates 



The X-ray ionization rate was investigated by Krolik &: Kallman (1983) and further developed 
by Glassgold et al. (1997). Glassgold et al. (2000) wrote the secondary electron contribution to the 
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X-ray ionization rate as follows: 



kT^ And? 



At 



J(t,xo), (11) 



where is the observed X-ray luminosity, and a{kTx) and T{kTx) are the total photoelectric 
absorption cross section and the optical depth at the energy E = kT^ respectively. The distance 
between the X-ray source and some point of the disk surface is denoted by "d" and the energy to 
make an ion pair is Ae. In above equation, the first factor between the square brackets corresponds 
to the primary ionization rate, assuming the same energy E = kTx for all primary electrons. The 
second factor kTx/ Ae then reads as the number of secondary electrons produced by a photoelectron 
with the energy kT^. The last factor J(r, xq) represents the attenuation of X-rays. Using the 
dimensionless energy parameter x = E/kT^, the attenuation factor J{t,xq) is written as 



/•oo 

j(^r,xo)= x-'^e-^-^('='^-)^""dx . (12) 

The optical depth T{kTx) measures how opaque the disk is toward the X-ray radiation from the 
star and is written as 

T{kT^) = Nh (T{kn) , (13) 

where a {kTx) is 

We assumed that heavy elements are depleted onto grains and used ctq = 8.50 x 10~^^ cm^ and 
n = 2.810 (Glassgold et al. 1997). The surface number density Nh is measured along the radiation 
path from the X-ray source. Letting a' be the angle between the radiation path and the radial 
axis, the surface number density is 



Nh = 



N, 



sm Q 

f oo 
iz 



J^n{a^^^ (15) 



where N±_ is the vertical surface density and n(a, z) is the number density at the disk radius "a" 
and the height "z". Note that n(a, z) and a' depend on the disk model (see eq. (17) and (21) 
respectively) . 

Usually, the ionization rate by cosmic rays is estimated as CcR,o — 10~^^s~^ (Spitzer & Tomasko 
1968), but it's not generally known to within better than an order of magnitude (Glassgold et al. 
2000). In some cases, however, the cosmic-ray ionization rate is well constrained. As an example, 
van der Tak & van Dishoeck (2000) obtained the cosmic-ray ionization rate of 2.6 it 1.8 x lO^^^s"^ 
through H^^CO submillimeter emission lines from massive protostars. This value is in good agree- 
ment with Voyager/Pioneer data (4 x 10~^^s~^). Since the attenuation length for cosmic rays is 
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XCR ~ 96 g cm ^ (Umebayashi &; Nakano 1981), we follow Sano et al. (2000) and write the cosmic 
ray ionization rate as follows 

S(a)-x(a,.)lA (^g) 
XCR \) 

where x(a, z) is the vertical column density measured at some height z (x(a, z) = p{a, z')dz'). 

The radioactive elements in the disk is yet another source of ionization, though they usually 
have only a minor effect on the total ionization rate. Following Umebayashi & Nakano (1981), we 
use the ionization rate by radioactive elements of (ra = 6.9 x 10~^^s~^. 



CCR 



CcR,i 



exp 



X{a,z) 
XCR . 



+ exp 



3.2. Disk Model 

Chiang &; Goldreich (1997), following Kenyon &; Hartmann (1987), developed a self-consistent 
passive disk model. Their model has a two-layer structure - a high temperature surface layer and 
a lower temperature interior. Dust grains in the surface layer of the disk absorb the flux of ultra- 
violet (UV) photons from the central star. Half of the emission from the grains in the surface layer 
escapes into the space and the remaining half heats the disk interior up. The model assumes the 
vertical hydrostatic equilibrium so that the number density may be expressed as 

n{a,z) = no(a) e , (17) 
where the number density at the midplane hq is 

no(a) = S(a)/(V27r/i(a)V5) • (18) 

Here, fig = 3 x 10~^^ g is the mean molecular weight of the gas, T,{a) is the surface density and 
h{a) is the pressure scale height (Chiang et al. 2001): 

S(a) = Eo (^)"'^' (19) 



where Sq is the surface mass density at 1 AU, Tc = GM^/ig/kR^ is a virial temperature that 
measures the gravitational potential at the surface of the central star, and and ii* are the stellar 
mass and radius respectively. Note that the disk interior temperature Tj (the disk temperature 
below the pressure scale height "h") decreases with a disk radius, so h{a) increases slightly weaker 
than the power of three halves (see eq. (20)) and no (a) decreases slightly weaker than the power 
of three (see eq. (18)). The pressure scale height h{a) is the height measured from the midplane 
to the interior disk, not the surface layer of the disk. 



The ionization of this disk by cosmic rays and radioactive elements is straight-forward to 
calculate. The X-ray ionization requires that we model the geometry of the source of the X-rays 
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- the central young stellar object. X-ray emission is thought to arise from the reconnection of 

large-scale magnetized loops, that are possibly but not necessarily related to the magnetopause 
radius of the stellar magnetic field and the inner edge of the disk. Accordingly, wc place the X-ray 
source in this picture at a fiducial distance of 12Rq above the disk and 12i?0 away from the central 
star following Glassgold et al. (1997). 

The X-rays from this source graze the surface of the flaring disk and penetrate it to produce 
the X-ray induced ionization. It is convenient to relate this grazing angle, = a', to the grazing angle 
defined in the disk model of Chiang et al. (2001) = a (see eq. (5) in Chiang et al. 2001). Assuming 
that the average radiation from the central star originates at a distance i?*/2 from the disk plane 
on the stellar surface (Hartmann 2000), we write the grazing angle from the X-ray source (see eq. 



where H is the height of the surface layer of the disk and defined as H/h = 4 in Chiang et al. 
(2001). Fig. 1 is a schematic figure of this equation. The grazing angle from the X-ray source, a' 
is the angle between the disk midplane and the solid arrow. The first term of the righthand side 
of the equation, a is the angle between the disk surface and the dashed arrow of the UV emission 
from the star. The second term f3 shows the flaring angle of the disk and the third term 7 comes 
from the geometry of the X-ray source. Note that 7 = corresponds to an X-ray source that is at 
the stellar surface. 

We calculated the X-ray ionization rate (eq. (11)) by integrating the attenuation factor (eq. 
(12)) using MATLAB. Instead of an infinite upper limit of the energy, we used x = 100 for the 
upper limit of the integration. For the lower limit, we used x = 1 keY /kT^ {E = 1 keV), following 
Glassgold et al. (1997). The surface number density along the path of the radiation Nh found in 
the optical depth (eq. (13)) is calculated by using the grazing angle a' (eq. (21)) and the number 
density n{a, z) that is obtained from the disk model of Chiang et al. (2001) (see eq. (17) and 
(18)). Note that this number density n(a,z) depends on the disk interior temperature Tj and the 
pressure scale height of the disk h. Both of them are determined by solving the radiative equilibrium 
equations (3) and (4) in Chiang et al. (2001). 



(15)) as 



a' = a — l3 + j 




(21) 



4. 



Results 



In this section, we study the extent of dead zones for various disks. We show that X-rays and 
cosmic rays are comparable ionization sources and determine the size of dead zones from the total 
ionization rate. Throughout this section, we assume the metal poor disk (set /3 = /J^ in eq. (9)). 
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4.1. Standard Disk Model 

First, we show the results of a fiducial model which uses the standard disk model of Chiang 
et al. (2001), the critical magnetic Reynolds number ReM,crit = 1 with aturb = 1) 0.1, 0.01, and 
when X-rays exist, the X-ray luminosity of = 10^^ erg with the temperature of kT^ = 1 
kcV. Wc chose the standard X-ray luminosity of Lx = lO^^erg s^^ because typical young stars have 
X-ray luminosities Lx ~ 10^^ — lO^^'erg s^^. Though the standard disk of Chiang et al. (2001) 
has the surface mass density at 1 AU of Sq = 10'^ g cm^^, we don't call this minimum mass solar 
nebular model because these two models are qualitatively different from each other. With the 
standard disk model of Chiang et al. (2001), we calculated the surface number density (eq. (15)) 
and/or the surface mass density (x in eq. (16)) to estimate the X-ray and/or cosmic ray ionization 
rate at each radius and height of the 2D disk. All disk parameters are the same as in Table 1 of 
Chiang et al. (2001). Other parameters used in this subsection is summarized in Table 2 in this 
paper. We determined the critical height of the dead zone at each radius of the disk from eq. (7). 

The size of dead zones arising solely from X-ray irradiation are plotted in Fig. 2. The figure 
shows the vertical cross section of the standard disk. The region below the dashed, long-dashed 
and dot-dashed lines show the dead zones. As the magnetic field becomes weaker (the alpha value 
changes aturb = Ij O-lj 0.01), the dead zone becomes larger. The result follows from the fact that 
as the magnetic field becomes weak, the growth rate of the MRI instability is reduced with respect 
to the local damping rate (see eq. (4) and the text therein). 

The equivalent result for cosmic rays is shown in Fig. 3. In this case, we just replaced the 
ionization rate in electron fraction (eq. (9)) with the cosmic ray ionization rate calculated by eq. 
(16). Again, we find the smaller dead zone as the magnetic field becomes stronger. 

These two figures show that cosmic rays dominate 1 keV X-rays in ionizing a disk. The dead 
zone stretches out to ~ 5 — 7 AU in the case of illumination by X-rays and only to ~ 1 — 2 AU in 
cosmic ray case. 

The dividing line between X-ray or cosmic ray ionized regions may be found by equating their 
ionization rates. We call this resulting vertical scale height the Rontgen height. Note that due to 
the higher number density, the X-ray ionization rate decreases toward the disk midplane as does 
the electron fraction, while the cosmic ray ionization rate is about the same everywhere in the disk. 
X-ray ionization dominates cosmic ray ionization above the Rontgen height. 

Fig. 4 shows the dead zones estimated by the total ionization rate. Again, we took a surface 
mass density at 1 AU of Sq = 10^ g cm~^, the dead zone criterion i?eM < 1, an X-ray luminosity of 
Lx = 10^^ erg s"-*^, and an X-ray energy of kT^ = 1 keV. The Rontgen height suggests that X-rays 
dominate cosmic rays only in the surface layer — where the disk is optically thin. Comparing Fig. 
3 with 4, it is apparent that the size of dead zones are determined by cosmic rays for the standard 
disk case. 

There are two characteristic dead zone radii. One of them is the midplane radius where the 
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dead zone surface cuts the disk midplane. The other one is the fiducial radius where the dead zone 

surface crosses the pressure scale height of the interior disk. We take this latter case to be our 
"fiducial" value because planets are likely to form in the interior disk that has much higher density 
(the order of ~ 10^) and larger particles (0.01 — 1000//m instead of 0.01 — l/im) compared to the 
surface layers. 

4.2. Sensitivity of Dead Zone Radius to Model Parameters 

Next, we compare cases of different parameters (Soj/cT^,) with one another. This is done to 
reveal which parameter has the largest effects on the size of dead zone. We also present results for 
two cases; one for only X-ray induced ionization, and the other for ionization arising from the total 
ionization rate. Wc make this distinction in order that the contribution of X-rays can be clearly 
discerned (cf. discussion at the beginning of §5). Throughout this subsection, wc set ReM,crit = 1 
and aturb = 0.01 (see Table 2). These are the conditions suggested by, for example, Gammie (1996). 

Fig. 5 is the plot of the fiducial dead zone radius estimated by X-ray alone as a function of 
the X-ray energy {kT^ = 1, 2, 3, 5, and 10 keV). For the X-ray ionization calculation, we used 
the typical X-ray luminosity of = 10^^ erg s~^. We include the minimum mass solar nebula 
(So = 10^ g cm~^) as well as a very massive disk of Sq = 10^ g cm~^ suggested by Murray et al. 
(1998) for Jupiter to migrate from 5 AU to < 1 AU. It is apparent that the sizes of dead zones 
decrease as the X-ray energy increases and/or the disk surface mass density decreases. 

Note that the figure depends on the choice of the critical magnetic Reynolds number ReM,crit, 
the alpha parameter aturb, and the X-ray luminosity Lx- For example, if we increase the X- 
ray luminosity by two orders of magnitude (Lx = lO^^erg s^^) in Fig. 5, then the resulting 
curves are pushed down by the factor of 10. This is because of the definition of the magnetic 
Reynolds number iigM = V aturftC [cgh/ {234y/ l3nT)]. Here, the bracket term is determined by a 
disk model (i.e. independent of the X-ray energy or luminosity) so that the magnetic Reynolds 
number is changed depending on at^rb a^^d For X-rays, the ionization rate is written as 
(^x = Lx [(j{kTx)/{4:TTd?Ae) J(r, zq)] where the luminosity changes the result by two orders of 
magnitude in our case (we used L^ = lO^^erg s~^ and lO^^erg s~^). 

The results for the dead zone radius for a disk undergoing the total ionization rate is shown in 
Fig. 6. It shows that the radii of dead zones are almost constant for X-ray energy of kTx = 1 — 3 keV 
and coincide with the results of Fig. 5 in the higher X-ray energy region. 

These figures also show the obvious point that the size of dead zones is not afi^cctcd by low 
energy X-rays (kTx ~ 1 — 3 keV) and that the dead zone gets smaller as the surface mass density 
decreases. The decrease of the surface mass density (from upper lines to lower ones) leads to a 
smaller optical depth r and a larger ionization rate (, so that the dead zone becomes smaller. 

We summarize the effect of changing parameters. If the X-ray luminosity is larger (lO'^^erg s~^ 
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rather than 10^^ erg s~^), X-rays penetrate deeper in the disk and therefore the dead zone becomes 
smaller. If the critical magnetic Reynolds number is larger (100 rather than 1), a very small 
magnetic field diffusivity can destroy the turbulence - the viscosity becomes negligible and therefore 
the dead zone tends to be larger. If the alpha parameter is smaller (0.01 rather than 1), then the 
magnetic field is weaker, turbulent viscosity is smaller and the dead zone becomes larger. In short, 
to have a smaller dead zone, a disk has to have a smaller surface mass density Sq, a larger X-ray 
luminosity and a larger aturh value as well as a smaller critical Reynolds number ReM,crit- 

To draw the reader's attention to our major result, we present Fig. 7 and 8 and study the 
effect of the surface mass density at 1 AU Sq. Fig. 7 uses the same data sets as Fig. 5, but 

here we plot the fiducial dead zone radius as a function of the surface mass density at 1 AU 
(So = 10^, 5 X 10^, 10^, 6 X 10^, and 10^ g cm~^) instead of the X-ray energy. It is apparent that 
if cosmic rays penetrate disks without being swept away, they dominate X-rays in ionizing a disk 
except for kT^ ~ 5 — 10 keV. 

Fig. 8 shows the equivalent result for the total ionization rate and uses the same data sets 
as Fig. 6. The X-rays have a rather small effect on the size of dead zones even when they have 
high energy {kT^ ~ 5 — 10 keV). The most striking result is that dead zone radii vary by at most 
an order of magnitude, from 1 — 10 AU for a range of two orders of magnitude in disk surface 
density and an order of magnitude in the X-ray energy. The robustness of this result indicates that 
the region of Jovian-mass planet formation is rather similar for a wide variety of protostellar disk 
systems. 

Wc can sec that the differences between lines in Fig. 7, 8 arc smaller than those in Fig. 5, 6. It 
appears that the size of the dead zone is most sensitive to the surface mass density and, to a lesser 
extent, the X-ray energy. These figures also show that the dead zone produced by cosmic rays is 
usually smaller than that by X-rays. X-rays dominate cosmic rays in ionizing the circumstellar disk 
only when the X-ray source has a high energy {kT^ ~ 5 — 10 keV). 

Table 3 shows the size range of dead zones for four extreme cases - {ReM,crit, otturb) = (!> 1)) 
(1, 0.01), (100, 1), and (100, 0.01). The minimum and maximum dead zone radii correspond to the 
minimum and maximum surface mass density at 1 AU. It is apparent that the dead zone is sensitive 
not only to the surface mass density or X-ray energy, but also to the value of the magnetic Reynolds 
number. The important point of this table is that the dead zone radii for total ionization rate are 
about the same as the region of terrestrial planets in our solar system except the one extreme case. 
For this extreme case of {ReM,critiOtturh) = (100,0.01), which is equivalent to R'^m crit ~ ^^^^ or 
{ReM,crit, oiturb) = (1; 10~^), the dead zone gets pushed out to 5.8 AU. 

Finally, we note that the data presented in Figures 5, 7 and 8 corresponding to the variation 
of dead zone radii can be well fit by approximate linear curves. For completeness, we present these 
fits in Appendix A - Table 5, 6, 7, and 8. 
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5. Application to the case of A A Tau 

We now specialize our model to an observed protostellar disk - AA Tau. We consider this 
source because of the excellent fit of the Chiang et al. model to the data from this system. The 
stellar and disk parameters that we used are listed in Table 4 and are taken from Chiang et al. (2001) 
other than X-ray parameters. The X-ray himinosity Lx and the X-ray energy kTx are taken from 
Neuhauser et al. (1995). For the dead zone criterion, we used ReM,crit = 1 with aturb = 1) 0.1, 0.01. 
These parameters are summarized in Table 2. 

Fig. 9 is the estimated dead zone for AA Tau. It shows that X-ray ionization dominates cosmic 
ray ionization only in the surface layer where the optical depth is very small (see Rontgen height). 
The dead zone becomes smaller as the magnetic field gets stronger because of the stronger viscous 
torque that it exerts. In all cases the dead zones by cosmic rays are smaller than those produced 
by X-rays. 

If the disk of AA Tau is metal dominant, the dead zone gets smaller. Fig. 10 shows the dead 
zones estimated by total ionization rate for both metal poor and metal dominant cases. Since 
metal ions recombine with electrons slowly, the electron fraction for metal dominant case is larger 
(compare Pd with f3r) - ReM tends to be larger - and hence the dead zone is smaller. Note that 
the recombination coefficient of metals is ~ 5 orders of magnitude smaller than that of molecular 
ions. The "real" dead zone would be calculated by solving eq. (8) directly and would be located 
somewhere between these two extremes. 

6. Discussion and Conclusions 

We calculated the dead zones for a variety of disk models, as well as the specific case of a 
Class II source, AA Tau. This kind of calculation has been done by several authors (e.g., Gammie 
1996; Igea & Glassgold 1999; Sano et al. 2000; Promang et al. 2002), but all of them use either the 

minimum mass solar nebula model developed by Hayashi et al. (1985) (e.g., Gammie 1996; Igea 
& Glassgold 1999; Sano et al. 2000) or the a disk model developed by Shakura & Sunyaev (1973) 
(Fromang et al. 2002). All our disk models are obtained by using the self-consistent passive disk 
model of Chiang et al. (2001). To calculate X-ray ionization rates, we followed the method used by 
Glassgold et al. (2000), while for cosmic rays, we used the method adopted by Sano et al. (2000). 
When the ionization rate becomes sufficiently low so that the magnetic Reynolds number is less 
than some critical value (1 or 100), we determined the size of the dead zone. We compared the 
result by X-rays with cosmic rays and also obtained the size of the dead zone determined by the 
total ionization rate. 

The fact that our dead zones extend as far as the region of terrestrial planet formation in our 
solar system is very important. In the accretion picture, considerable gas accretion must occur 
onto a sufficiently massive protoplanetary core. This must inevitably take place in the turbulent 
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region of the disk beyond the dead zone. The major imphcation of our work is that the division 

between Jovian and sub- Jovian or terrestrial planets may be due to the presence of a dead zone. 
Wc explore the significance of these results for predicting planetary masses in a forthcoming paper. 
Our basic results are as follows: 

1. Our major finding is that the typical dead zone encompasses, in physical scale, the terrestrial 
planets in our solar system. 

A typical dead zone size estimated from the total ionization rate is ~ 0.24 — 2.7 AU for 
the fiducial disk. For an extreme case of {ReM,crit-,Oiturb) = (100,0.01) {R'^m crit ~ 1000 or 
{ReM,crit, Oiturb) = (1) 10~^) equivalently) , the dead zone extends out to ~ 5.8 AU. These re- 
sults suggest that the observed exosolar Jovian-mass planets must have migrated to their observed 
positions from points of origin farther out in disk radius, beyond these dead zone radii. 

2. X-rays dominate cosmic rays in ionizing a disk only in the surface of the disk for most cases. 

Glassgold et al. (2000) noted that X-ray ionization dominates cosmic ray ionization out to 
~ 1000 AU. We found that the disks are usually too optically thick for this to be generally true. 
X-rays only dominate cosmic rays in the surface layer of the disk. Cosmic rays determine the size 
of the dead zone under typical conditions (compare, for example, Fig. 3 with Fig. 4). 

3. Sufficiently high energy X-rays could dominate cosmic rays in ionizing a disk. 

The X-rays could dominate the cosmic rays if the X-ray energy is high kT^ = 5 — lOkeV 
(compare Fig. 7 with 8). This is much higher energy for most observed sources however. 

4. The size of a dead zone is sensitive to the disk surface mass density Eq and the X-ray energy 

We found that the size of the dead zone is most sensitive to the surface mass density, and to a 
lesser extent, the X-ray energy (see Fig. 5 & 6 and Fig. 7 & 8). The X-ray luminosity, the critical 
magnetic Reynolds number, the alpha viscous parameter however, have a significant effect on the 
size of a dead zone. 

5. There is a power law relation between the size of a dead zone estimated by X-rays alone 
and So or kT^ (see Appendix). 

The dead zone radius has power-law relations with both the surface mass density and the X-ray 
energy. The dead zone radii estimated by X-rays obey oc Sg^^ — Eg^^ and oc kT~'^''^^ — kT~^'^ , 
while those estimated by cosmic rays obey ad oc - The dead zone radii determined by 

the total ionization rate has the relation of aa oc Sq^® — Eg^^. 

We thank Norm Murray, Joe Weingartuer and Steve Balbus for stimulating conversations on 
these topics. We also thank an anonymous referee for useful review of the manuscript. SM is 
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Table 4: Input parameters for AA Tau 



Symbol 


Meaning 


Standard value 


M, 


stellar mass 


O.67M0 


R* 


stellar radius 


2.1Rq 




stellar effective temperature 


4000 K 


Lx 


X-ray luminosity 


0.439 X 10=^° erg 


kTx 


X-ray energy 


1.21 keV 


T-iiion 

sub 


iron sublimation temperature 


2000 K 


rpsil 

sub 


silicate sublimation temperature 


1500 K 


rpice 

sub 


ice H2O sublimation temperature 


150 K 


So 


surface mass density at 1 AU 


1.5 X 10^ g cm~^ 


P 


— dlog S/dloga 


1.5 


Qi 


inner disk radius 


2 i?* 


tto 


outer disk radius 


8600 i?* = 250 AU 


H/h 


visible photosphere height / gas scale height 


3.8 


Qi 


—dlogN/dlogr in the interior disk 


3.5 


Qs 


—d log N/ d log r in the surface disks 


3.5 




maximum grain radius in the interior disk 


1000 nm 




maximum grain radius in the surface disks 


1 fim 
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Appendix A 

X-rays and Cosmic rays 

We supply tables of X-ray and cosmic ray power-law relations separately for interested readers. 
These approximate relations are all obtained by using the least square fitting. 

Table 5 shows the fitting parameters for the dead zone radius by X-rays as a function of X-ray 

energy kT^. For the same choice of the critical magnetic Reynolds number ReM,crit and the alpha 
parameter aturh-, it is apparent that the curves show roughly the same power. One of the examples 
is shown in Fig. 5. 

Tabic 6 shows the fitting parameters for the dead zone radius by X-rays as a function of the 
surface mass density at 1 AU, Sq. Again, we can see the similar power for the same choice of the 
critical magnetic Reynolds number ReM,crit and the alpha parameter aturb- One of the examples 
is shown in Fig. 7. 

Table 7 shows the fitting parameters for the dead zone radius by cosmic rays as a function of 
the surface mass density at 1 AU, Sq. One of the examples is shown in Fig. 7 (the gray solid line 
with crosses). 

Table 8 shows the fitting parameters for the dead zone radius estimated by the total ionization 
rate as a function of the surface mass density at 1 AU, Sq. Again, we can see the similar power for 
the same choice of the critical magnetic Reynolds number ReM,crit and the alpha parameter aturb- 
One of the examples is shown in Fig. 8. In particular, kT^ = 1 — 3 keV give a very similar power 
law relation. 
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Table 5: The sizes of dead zones as a function of kT^: y = ax . The star has X-ray luminosity of 
10^^ erg s~^. 



{ReM,crit-i Oiturb) 


So[g cm-2] 


a 


h 




(1,1) 


10^ 


2.67 


-1.05 


0.692e 


- 4 




5 X 10^ 


5.97 


-0.908 


0.583e 


- 3 




10^ 


8.47 


-0.887 


0.140e 


- 2 




6 X 10"^ 


27.3 


-0.887 


0.140e 


-2 




10^ 


28.6 


-0.883 


0.164e 


-2 


(1,0.01) 


10^ 


3.73 


-0.776 


0.179e 


-2 




5 X 10^ 


8.58 


-0.761 


0.201e 


- 2 




10^ 


11.9 


-0.751 


0.213e 


- 2 




6 X 10^ 


29.7 


-0.772 


0.267e 


- 2 




10^ 


38.5 


-0.786 


0.238e 


-2 


(100,1) 


10^ 


5.86 


-0.609 


0.403e 


-2 




5 X 10^ 


12.8 


-0.619 


0.321e 


-2 




10^ 


18.4 


-0.640 


0.320e 


- 2 




6 X 10^ 


44.8 


-0.665 


0.370e 


- 2 




10^ 


57.9 


-0.676 


0.296e 


- 2 


(100,0.01) 


10^ 


9.59 


-0.458 


0.450e 


-2 




5 X 10^ 


20.8 


-0.483 


0.503e 


-2 




10^ 


29.4 


-0.508 


0.383e 


-2 




6 X 10^ 


69.7 


-0.529 


0.387e 


-2 




10^ 


89.5 


-0.544 


0.423e 


-2 
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Table 6: The sizes of dead zones as a function of Eq: y = ax". The star has the X-ray luminosity 
of lO^^erg s~^. 





/,-r,[k(-v] 


n 


b 


9 

Y" 


(1,1) 


1 


0.761e - 


1 


0.515 


0.315e 


-4 




2 


0.359e - 


1 


0.522 


0.443e 


- 3 




3 


0.225e - 


1 


0.530 


0.658e 


- 3 




5 


0.123e - 


1 


0.545 


0.220e 


- 2 




10 


0.454e - 


2 


0.589 


0.820e 


- 2 


(1,0.01) 


1 


0.118 




0.508 


0.483e 


-4 




2 


0.676e - 


1 


0.501 


0.977e 


- 4 




3 


0.488e - 


1 


0.500 


0.258e 


- 3 




5 


0.330e - 


1 


0.502 


0.477e 


- 3 




10 


0.213e - 


1 


0.501 


0.757e 


-3 


(100,1) 


1 


0.201 




0.498 


0.496e 


-4 




2 


0.131 




0.483 


0.178e 


-5 




3 


0.995e - 


1 


0.483 


0.1 78e 


- 5 




5 


0.775e - 


1 


0.476 


0.474e 


-4 




10 


0.652e - 


1 


0.460 


0.324e 


-4 


(100,0.01) 


1 


0.369 




0.483 


0.178e 


- 5 




2 


0.256 




0.471 


0.158e 


-4 




3 


0.212 




0.467 


0.387e 


-4 




5 


0.187 




0.455 


0.469e 


-4 




10 


0.171 




0.441 


0.723e 


-4 



Table 7: The sizes of dead zones produced by cosmic rays, as a function of Sq: y = ax". 



{ReM,critj (^turb) 


a 


b 


x' 


(1,1) 


0.234e - 1 


0.485 


0.978e 


-3 


(1,0.01) 


0.650e - 1 


0.420 


0.160e 


-2 


(100,1) 


0.236 


0.341 


0.209e 


- 2 


(100,0.01) 


0.830 


0.278 


0.411e 


-3 
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Table 8: The sizes of dead zones estimated from the total ionization rate as a function of Sq: 
y = ax''. The star has the X-ray luminosity of 10^^ erg s~^. 



{ReM,criti (^turb) 




a 


b 


x' 


(1,1) 


1 


0.234e - 


1 


0.485 


0.978e - 


3 




2 


0.234e - 


1 


0.485 


0.978e - 


3 




3 


0.217e - 


1 


0.492 


0.842e - 


3 




5 


0.125e - 


1 


0.539 


0.854e - 


3 




10 


0.454e - 


2 


0.589 
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Fig. 1. — A schematic figure of the cross section of the disk. The soUd arrow shows the X-rays 
from the magnetic field (the X-ray source is at (12i?o, I2-R0)) and the dashed arrow shows the 
hght from the central star. Note that there is a relation of a' = a — /? + 7 among angles (see eq. 
(21) and the text). 
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Fig. 2. — The dead zone estimated by X-rays for the standard disk model of Chiang et al. (2001). 
We used Lx = 10^^ erg s~^, kTx = 1 keV, ReM,crit = 1 and aturb = 1; 0-1) ^ 0.01. X-ray dead 
zone stretches out to 5 — 7 AU. 
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Fig. 3. — The dead zone estimated by cosmic rays for the standard disk model of Chiang et al. 
(2001). We used ReM,crit = 1 and amrb = li 0.1, & 0.01. Cosmic ray dead zone stretches out to 
1 - 2 AU. 
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Fig. 4. — The dead zone determined by the total ionization rate for the standard disk model 
of Chiang et al. (2001). We used Lx = 10^^ erg s~^, kT^ = 1 keV, ReM,crit = 1 and aturb = 
1, 0.1, & 0.01. Also shown is the Rontgen height. X-rays dominate cosmic rays in ionizing a disk 
only at the surface layer, so the dead zone estimated from the total ionization rate has the same 
size as that from cosmic rays. 
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Fig. 5. — The fiducial dead zone radius as a function of the X-ray energy for a source luminosity 
Lx = 10^^ erg s~^. We chose ReM,crit = 1 and aturb = 0-01 ^^nd estimated the dead zone sizes by 
X-rays alone. Higher X-ray energy and smaller surface mass density lead to the smaller dead zone. 
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Fig. 6. — The fiducial dead zone radius as a function of the X-ray energy for a source luminosity 
Lx = 10^^ erg s~^. We chose ReM,crit = 1 and aturb = 0-01 and estimated the dead zone sizes by 
total ionization rates. Cosmic rays dominate low energy X-rays in ionizing a disk. 
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Fig. 7. — The fiducial dead zone radius as a function of the surface mass density for an X-ray 
luminosity = 10^^ erg s~^. We chose ReM,crit = 1 and aturb = 0.01 and estimated the dead zone 
sizes by X-rays alone. Also shown is the dead zone radius estimated by cosmic rays. High energy 
X-rays dominate cosmic rays in ionizing a disk. 
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Fig. 8. — The fiducial dead zone radius as a function of the surface mass density for an X-ray 
luminosity = 10^^ erg s~^. We chose ReM,crit = 1 and aturb = 0.01 and estimated the dead zone 
sizes by total ionization rates. The dead zone radius changes about an order of magnitude for a 
range of two orders of magnitude in disk surface mass density. 
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Fig. 9. — The dead zone of A A Tau estimated from X-rays is compared with that predicted from 
cosmic rays. We used = 0.439 x 10^'' erg s^^, kT^ = 1.21 keV, ReM,crit = 1, and aturb = 
1, 0.1, and 0.01. Also shown is the Rontgen height. Cosmic rays dominate X-rays in ionization 
almost throughout the disk. 
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Fig. 10. — The dead zone of AA Tau estimated both for a metal poor and a metal dominant disk. 
We used = 0.439 x 10^° erg s'^ , kT^ = 1.21 keV, ReM,crit = 1, and aturh = 1, 0.1, and 0.01. 
The presence of metals makes the dead zone smaller. 



